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The electronic energy levels of one-dimensional aperiodic systems driven by a homogeneous electric 
field are studied by means of a phase space description based on the Wigner distribution function. 
The formulation provides physical insight into the quantum nature of the electronic states for the 
aperiodic systems generated by the Fibonacci and Thue-Morse sequences. The nonclassical param- 
eter for electronic states is studied as a function of the magnitude of homogeneous electric field to 
achieve the main result of this work which is to prove that the nonclassical properties of the elec- 
tronic states in the aperiodic systems determine the transition probability between electronic states 
in the region of anticrossings. The localisation properties of electronic states and the uncertainty 
product of momentum and position variables are also calculated as functions of the electric field[^ 
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I. INTRODUCTION 

Analysis of the electronic states of artificial structures 
(e.g. superlattices, quantum wires, rings or dots) plays 
a central role in modern condensed matter physics be- 
cause they determine many useful properties indispens- 
able in industrial applications of atomic- or nano-scale 
devices. The experimental realisation of that kind of 
structures allows to investigate the influence of quan- 
tum effects on mechanical, optical or transport prop- 
erties. The atomic cluster structures such as quantum 
corrals or chains can be fabricated using the low temper- 
ature scanning tunnelling microscope (STM) manipula- 
tion of individual atoms on the conducting substrat^^^ 
or mechanically controllable break junction (MCBJjP'^. 
The latter method allows to fabricate only very short 
mono or mixed atomic chains that consist of 4 or 5 metal 
atomd^Zl on the contrary, the tip of STM can be used 
to build longer and much more complex structures with 
different shapes. Additionally, long atomic chains can be 
reconstructed on a flat surface by self-assembly. Exper- 
imental results for self-assembled gold atomic chains on 
silicon surfaces conflrmed the existence of long and sta- 
ble chain#12l this case the chains are extended over 
hundreds of nanometres and can break due to atomic 
defects or because of the intentional removal of single 
atom or group of atoms from the perfect chain. It means 
that experimental techniques open up the possibility to 
fabricate novel metal nanostructures created by the in- 
tentional atomic rearrangement. On the other hand, the 
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external fields interacting with electronic states of the 
systems can modify their electronic spectrum and new 
properties of the system are observed. Especially the in- 
teraction with the electronic states in the system with 
broken translational symmetry (disordered or aperiodic 
systems) seems to be interesting because the quantum 
interference effects determine the electronic properties 
of the systems at low temperatures e.g.liSHlll, One of 
the simplest but nontrivial examples of the interactions 
is the effect of the homogeneous electric field on elec- 
tronic states in crystals where the Wannier-Stark quanti- 
sation was predictediSJMI^ g^^^ confirmed experimentally 
in artificial semiconductor and optical superlatticesl^HlSl, 
Discussion of the effect in one-dimensional periodic sys- 
tems has been carried out for many years^SHSH jg gt 
the subject of current experimental as well as theoretical 
research^^^^^'. 

In this paper, we apply the phase space approach 
based on the non-classical distribution functionp2H36| to 
the problem of electronic states in isolated and finite 
one-dimensional aperiodic systems generated by the Fi- 
bonacci and Thue-Morse sequenceP^HSS] jjj the presence 
of an external homogeneous electric field. This ap- 
proach is widely used in the study of quantum trans- 
port phenomen£p2lHll^ ^j^j^t less exploited in the descrip- 
tion of electronic states of nanosystem^^. The phase 
space method allows one to carefully investigate the in- 
creasing role of the quantum effects resulting from inter- 
play between the electric field and aperiodic ordering in 
the one-dimensional systems and leading to the greater 
importance of the nonclassical properties of electronic 
states. We compute the Wigner distribution functions 
for some electronic states in the aperiodic systems and 
analyse their localisation properties in the phase space as 
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a function of electric field. The analysis is based o n the 
inverse participation ratio (IPR) in the phase spac^^^. 
This approach allows us to show that the electric field 
increases the nonclassical properties of electronic states 
close to the transition between them. 

The presented analysis can be simply extended to the 
description of cold atoms in the optical latticd^^^^ or 
layered systems such as aperiodic superlattices consisting 
of two different materials that are arranged according to 
appropriate aperiodic sequence, e.g.l^^l^^ and references 
therein. 

The rest of the paper is organised as follows. In Sec. 
we present the theoretical model of pseudoatomic chain. 
In Sec. |III| we apply the formalism of non-classical distri- 
bution functions to the quantitative analysis of electronic 
states in the aperiodic systems in the external electric 
field. In concluding remarks we summarise presented re- 
sults. 
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FIG. 1. Energies of states in the conduction band for N=100 
wells as a function of the electric field for (a) periodic, (b) 
Fibonacci and (c) Thue-Morse sequence. 



II. THEORETICAL MODEL AND METHODS 

We consider the aperiodic chain of metallic potential 
wells modelling atoms which can be described by the one- 
particle Hamiltonian in the form 



Ho 



2m 



N 



(1) 



where m is the effective mass of electron, and N is the 
number of wells. The positions of wells, Xi, are dis- 
tributed according to the binary Fibonacci and Thue- 
Morse sequences. We generate these sequences over set 
{0, 1} using the following inflation rules^ 

a. Fibonacci sequence: — > 01, and 1 — > 0, 

b. Thue-Morse sequence: — > 01, and 1 — > 10. 

In our notation zero corresponds to an empty site in the 
simple crystal lattice, and one corresponds to a site in 
the lattice occupied by the potential well v{x — Xi). 

The potential term in the Hamiltonian ([T]) is repre- 
sented by the superposition of the potential wells v{x — 
Xi). We assume that each well in the chain is given by 
the Shaw pseudopotential modified by screening, namely 



\ x—X \ ' I * I <^ 



CXp {—^Xc) 



\X - XA < Xr. 



(2) 



A quantity Xc is the cut-off parameter and it is chosen to 
have the value for which the pseudopotential reproduces 
the ionisation energy of Cu; /i is the screening factor in 
the Thomas-Fermi approximation. We assume that each 
well gives only one state to the conduction band and 
therefore the Fermi level, Ep, is defined in the middle of 
the conduction band. Finally a constant electric field £ is 



applied along the wire and the total Hamiltonian of the 
system under the electrostatic perturbation has a form 



H ^Hq + e£x, 



(3) 



where — e is the electron charge. 

The energy spectrum of the finite system which is de- 
scribed by the Hamiltonian ([T]) forms the energy bands. 
The structure of these energy bands strongly depends 
on the arrangement of potential wells. When a peri- 
odic system is subjected to the electric field, the eigen- 
states of the Hamiltonian ^ form so-called homogeneous 
Wannier-Stark ladders ^^ — '. This situation is presented 
in Fig. [T] a. In the aperiodic structures generated by 
the Fibonacci and Thue-Morse order, the eigenstates of 
the Hamiltonian ([s]) form more complex energy spectra 
of the conduction band as it is shown in Figs, [ij b and 
[T] c. In these cases, the electronic states tend to group 
into subbands. The energy widths of these subbands are 
much smaller than the width of the conduction band for 
the periodic system and strongly depend on the values 
of electric field and number of wells in the systems. As 
a result, the aperiodic systems form the inhomogeneous 
Wannier-Stark ladders. From these results we conclude 
that the effect of inhomogeneous Wannier-Stark ladders 
originates from the interplay between aperiodic order and 
the electric field on the electronic states. A closer inspec- 
tion of the inhomogeneous Wannier-Stark ladders reveals 
the occurrence of changes in energy levels inclinations in 
the electric field. These changes result from the repul- 
sion between neighbouring electronic states, and the an- 
ticrossingJ^ are observed between them for some values 
of the electric field as it is shown in Fig. [2] 

These non-trivial properties of the inhomogeneous 
Wannier-Stark ladders suggest that the correlated po- 
sition and momentum behaviour of electronic states play 
an important role in the region of anticrossings. There- 
fore the nature of these electronic states is analysed by 
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FIG. 2. Energies calculated for (a), (b) the Fibonacci se- 
quence and (c), (d) Thue-Morse sequence, (a) and (c) show 
states at the Fermi level (Ef) and directly below and above 
the Fermi level {Ep-i and -Ef+i, respectively), (b) and (d) 
show states at the bottom of the conduction band. Arrows 
point to the anticrossings. 

the phase-space methods. Here, we restrict ourselves to 
the detailed analysis of the states in the vicinity of the 
band bottom, but we also present briefly some results for 
the states in the middle of the conduction band. 

Within the phase space approach, the electronic states 
may be represented by the Wigner distribution function. 
For an electron in the pure state the Wigner distribution 
function has the fornPHEH 

fn{x,k) - J dx'{x^^x'\lbn){^n\x+^x')e''''''. (4) 

It should be noted that the Wigner distribution function 
is a bilinear combination of the electron wavefunction and 
therefore contains interference information. By means 
of the Wigner distribution function we can evaluate the 
expectation value of any Hcrmitian quantu m- me chanical 
operator A in the state n using the formulcP2125l 

(A)„ = J dxdkA{x,k)Uix,k), (5) 

where A{x, k) is the Wigner representation of quantum- 
mechanical operator A given by 

A{x, k) = J dx'{x ~ ^x'\A\x + ia;')e**=^'. (6) 

Because the Wigner distribution function can take nega- 
tive values in some subregions of the phase space it can- 
not be interpreted as the classical distribution function 



in the phase space. The negative part of the Wigner 
distribution function is responsible for quantum correla- 
tions between spatially separated pieces of the electronic 
stat^^. It stems from the fact that the information from 
the off-diagonal terms in Q represented by a;' 7^ vari- 
able is transferred to the Wigner distribution function 
via the momentum k. These properties characteristic for 
the Wigner distribution function may be utilised as an 
indicator of nonclassicality of electronic state^^^*^. 



III. RESULTS AND DISCUSSION 

In the present considerations we assume that the aperi- 
odic systems generated by the Fibonacci and Thue-Morse 
sequences are limited by the infinite wells. The energy 
spectra were calculated for weak electric fields using the 
parameters of the conduction electron in copper. In our 
calculations we consider chains of different lengths (e.g. 
30, 50 or 100 potential wells). Here we present the re- 
sults for 100 wells containing all characteristic properties 
of the obtained results. All values are given in atomic 
units ?i = |e| = — 1. 

In the first step we determine the Wigner distribution 
functions for some of pure states and investigate their 
changes due to the electric field. Fig. [3] shows the Wigner 
distribution functions for the two lowest electronic states 
in the Fibonacci chain in the vicinity of the first anti- 
crossing that is pointed to by the arrow in Fig. [2j b. 
Initially, both electronic states which are represented by 
the Wigner distribution functions fo{x,k) and fi{x,k) 
occupy different regions of the phase space. Increasing 
the electric field shifts the states in opposite directions 
so that the distance between them decreases. Finally, in 
the region of anticrossing, the Wigner distribution func- 
tions occupy the same region of the phase space, and it 
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FIG. 4. The transition probability between the states n and 
n + 1 at the bottom of the conduction band in the system 
based on the Fibonacci sequence; (a) n = 0, (b) n = 1. 
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FIG. 5. The Wigner function for the system based on the 
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leads to overlapping of these functions as shown in Figs 
[Sj b and|3] e. In this case we may express the electronic 
states by the linear superposition of individual states: 
W = aolV'o) + ailV^i), where |aoP + jaip = 1. Using 
the definition given by Eq. Q we obtain expression for 
the Wigner distribution function of the electronic state 
in the form 

f{x,k) = \ao\^fo{x,k) + \ai\^h{x,k) (7) 
+ 2K|aoaI J da:'(x - |^|V'i)(^oN + |^)e^'="'|. 

The first two terms correspond to individual electronic 
states, and the last term represents the mixture part 
which results from the bilinearity of the Wigner distri- 
bution function. Th is typ e of term is often termed the 
quantum interferencePsSl. 

Quantitative analysis of the transition between elec- 
tronic states Itpn) and \ipm) may be based on the prod- 
uct of their Wigner distribution functions fn(x,k) and 
fmix,k) integrated over the phase space^^"^, namely 

Pnm = '^TT J dxdk fnix,k) f^{x,k). (8) 

The results of calculations of the overlapping integral 
given by Eq.(|8]) for the three lowest states in the Fi- 
bonacci chain are shown in Fig. |4j The maximum values 
of the overlapping integrals in both cases correspond to 
the electric field values at which the anticrossings be- 
tween the appropriate states are observed. After a fur- 
ther increase in the electric field the quantum interference 
term disappears because the distance between states be- 
comes larger and finally both states are well separated in 
the phase space which means that the overlapping inte- 
grals tend to zero. A similar situation is observed for the 
electronic states in the Thue-Morse chain, for example in 
case of the Wigner functions shown in Fig. [5] 



A more complex situation is presented in Fig. |6]where 
the evolution of the Wigner distribution functions for 
three quantum states in the vicinity of the middle of the 
conduction band is shown as a function of the electric 
field. In this case the energy spectrum (see Fig. [2]) ex- 
hibits two anticrossings located close to each other. Each 
of these anticrossings mixes only two neighbouring elec- 
tronic states, therefore this process can be explained by 
the previous analysis. As it is presented in Figs [3] [5] [6] 
quantum states in the aperiodic chains under the fixed 
boundary conditions localise in limited areas of the phase 
space. The extents of these areas depend on the assumed 
distribution of the potential wells (based on Fibonacci or 
Thue-Morse sequence) and the number of quantum state. 

One of the possibilities of measuring the degree of lo- 
calisation is to u se the IPR parameter calculated in phase 
space, defined hj^^ 

lPRn = ^ J dxdk[g,,{x,k;l/2)]\ (9) 

where gn{x,k; 1/2) is the Husimi function. The Husimi 
function is an example of the non-negative non-classical 
distribution functions which can be obtained by the con- 
volution of the Wigner distribution function and a win- 
dow functiorP2l 

gn{x,k;A^k) ^ J dx'dk' W{x-x' ,k-k';A^k)fn{x\k'), 

(10) 

where W{x — x' ,k — k'; A^k) is a window function with 
resolution A^k- In particular, the Husimi function is ob- 
tained if we choose the window function as a Gaussian 
function with the resolution corresponding to the min- 
imum resulting from the uncertainty principle {A^k = 
1/2). 

IPR can be used to study the influence of the electric 
field on localisation, as shown in Fig. [Tj For the ground 
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FIG. 6. The Wigner function for the system based on the 
Fibonacci sequence, (a), (b), (c) - Ep-i, (d), (e), (f) - Ef, 
(g)-(i) - Ef+1. Electric field, from left to right, 1.8, 1.85, 1.9 
X 10"*' a.u. 
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FIG. 7. Inverse Participation Ratio ( IPRn) for (a), (b) the 
Fibonacci sequence and (c), (d) Thue-Morse sequence. 
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FIG. 8. The average position {x)n and squared momentum 
{k^)n for the Fibonacci sequence. 



state in the Fibonacci chain (see Fig. [7j b) it was found 
that after IPR reaches the value of 0.13, further increase 
in the electric field to 2 x 10~^ a.u. does not change the 
localisation. In a similar way, for the first excited state 
IPR stabilises initially at 0.08 for the fields larger than 
0.7 X 10~^ a.u. Then, when the electric field is between 
2 and 2.5 x 10~^ a.u., we observe that the energy levels 
values plotted versus the electric field change their gradi- 
ent (see Fig. [2] b). This behaviour is connected with the 
anticrossing taking place between the mentioned values 
of the electric field. In the phase space, the shape of the 
states changes, as well as the occupied area. Both states 
are shifted and the degree of localisation changes signifi- 
cantly. As a consequence, the first excited state localises 
in the region of the phase space previously occupied by 
the ground state, and vice versa. The expectation val- 
ues of the position and squared momentum presented on 
Figs |8] and [9] were calculated for these states. The results 
confirm that the electric field shifts the states only at the 
points of anticrossings. 

In case of the Thue-Morse sequence the observed be- 
haviour is slightly different. When the electric field is 
applied, localisation of the two lowest states starts to 
increase very quickly and the states are shifted. Then, 
the degree of localisation falls down and finally settles 
at a constant level. Also the expectation values of posi- 
tion and momentum do not undergo any further changes 
when the electric fields are increased. Similar character- 
istics of the ground state and the first excited state mean 
that increasing the electric field separates energies of the 
states, as shown in Fig. [2j d. 

Figs[7| a, c present changes in localisation of the states 
from the middle of the band. The changes are very sharp 
and are observed for the electric field values at which the 
plot of energy of the state against the electric field al- 
ters its gradient. Modification of the phase space areas 
occupied by the states observed at the anticrossings, to- 
gether with the change in IPR values, are caused by the 
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FIG. 9. The average position {a;)„ and squared momentum 
{k^)n for the Thue-Morse sequence. 



decreasing of the energy between the states. For that 
reason we expect the quantum effects to have a greater 
impact on the analysed states, notably the quantum in- 
terference between them. Appearance of that kind of 
effects should be accompanied by an increase in the neg- 
ative part of the Wigner function which is responsible for 
the quantum interference. 

To investigate this phenomenon we split up the Wigner 
distribution function into two parts, namely 



fnix, k) = f+{x, k) + /„ {x, fc). 



(11) 



where f^{x,k) and f~{x,k) correspond to the positive 
part and the negative part of the Wigner distribution 
function, respectively. 

The measure of the non-classical n ature of an electronic 
state \ipn) is defined by the formulgP^EH 



1 



1+ 



(12) 



where X+ and I~ are the moduli of the integrals of the 
positive part and the negative part of the Wigner dis- 
tribution function, respectively. Influence of the electric 
field on the f„ parameter is shown in Fig. [lO] 
In the Fibonacci chain the non-classical character consid- 
erably decreases (to about 0.52) immediately after turn- 
ing the electric field on. Then its value does not change 
until the anticrossing takes place in the energy spectrum, 
where the f„ parameter rises as a result of the increas- 
ing Wigner function negative part for both states. For a 
further increase in electric field the non-classical nature 
decreases. 

For the ground state in the Thue-Morse chain a mono- 
tonic decrease in the parameter Vn is observed. It can be 
explained on the basis of large separation in the phase 
space between the ground state and the excited states. 
States placed in the middle of the band also exhibit jumps 
of the non-classical behaviour at the electric fields at 



0.6 



0.4 



i 0.2 - 



0.6 - 



0.4 - 



(a) 



(b) 




(d) 



3 1 
Electric field [10"'' a.u.] 



FIG. 10. The nonclassicality parameter calculated for the 
ground state (n — 0) in case of (b) the Fibonacci sequence and 
(d) Thue-Morse sequence. The nonclassicality parameter v„ 
calculated for the states corresponding with the Fermi level 
in case of (a) the Fibonacci sequence and (c) Thue-lVIorse 
sequence. 



which the anticrossings are observed. Here, however, the 
widths of anticrossings are very small, and the numerical 
accuracy does not allow to illustrate the peaks of non- 
classicality parameter, but only to show the values below 
and above anticrossings. 

Additionally we calculate the standard deviations of 
position and momentum variables that are defined by ex- 
pressions: al{x) = {x'^)n-{x)l andcr2(fc) = {k'^)„-{k)l, 
respectively, applying Eq. ([5| to find the relevant quan- 
tities from the Wigner function fn{x,k). The acquain- 
tance of these deviations allows one to calculate the un- 
certainty product an{x)(Jn{k) as a function of the electric 
field. Fig. 11 shows results for the uncertainty product 
for the two lowest states in both aperiodic chains. We 
may see that the electric field produces rapid changes 
in the uncertainty product around the anticrossings for 
both states. This behaviour of the uncertainty product 
can be explained by the highly nonclassical features of the 
Wigner distribution function corresponding to the quan- 
tum state in the region of anticrossing where the quan- 
tum interference phenomena play important role. On the 
other hand the uncertainty product allows one to distin- 
guish between pseudoclassical and nonclassical states. It 
results from the fact that the uncertainty product is equal 
1/2 for the coherent states which are recognised as the 
states closest to classical ones. Hence the nonclassical 
states may be characterised by the uncertainty product 
greater than 1/2. As it presented in Fig. 11 the uncer- 
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FIG. 11. The product of the position and momentum vari- 
ances for the Fibonacci sequence (a) and the Tue-Morse se- 
quence (b). 

tainty product is usually greater than one which means 
that the discussed states are nonclassical. 



IV. CONCLUDING REMARKS 



We have applied the phase space method based on the 
Wigner distribution function to analyse the region of an- 
ticrossings for the lowest electronic states in both ape- 
riodic systems. The discussion of the electronic states 
in the middle of conduction band is also included. This 
formulation allows to investigate the nonclassical prop- 
erties of the electronic states and their influence on the 
transitions between them. We have shown that the non- 
classical properties of electronic states of the aperiodic 
systems under the homogeneous electric field are most 
profound in the regions of anticrossings. The increase in 
the nonclassical parameter in these regions is a conse- 
quence of the increasing role of the negative part of the 
Wigner distribution function and it is correlated with the 
changes in the localisation properties of electronic states 
and the dynamical variables. Finally, we have used the 
uncertainty product of momentum and position variables 
as a simple measure of nonclassicality of the electronic 
state. We have found that this product as a function of 
the electric field is also correlated with the nonclassical 
properties of electronic states in the aperiodic systems. 
In the paper we have paid attention to the aperiodic sys- 
tems in the limit of low homogeneous electric field. The 
limit of strong electric field and dynamical aspect of the 
problem will be included in forthcoming publications. 



We have studied the influence of the homogeneous 
electric field on the energy spectrum of the finite one- 
dimensional aperiodic systems generated by the Fi- 
bonacci and Thue-Morse sequences. We have shown that 
the homogeneous electric field modifies the energy spec- 
trum of aperiodic systems and leads to the formation of 
the subbands. The appearance of subbands is a conse- 
quence of the attraction between energy levels for some 
value of the electric field. The anticrossings are observed 
for these values of the electric field. 
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